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A note on functional limit theorems for compound Cox processes* * * § 

V. Yu. Korolev! A. V. Chertokf A. Yu. Korchagin^ E. V. Kossova^A. I. Zeifmaih 


Abstract: An improved version of the functional limit theorem is proved establishing weak convergence of random 
walks generated by compound doubly stochastic Poisson processes (compound Cox processes) to Levy processes in the 
Skorokhod space under more realistic moment conditions. As corollaries, theorems are proved on convergence of random 
walks with jumps having finite variances to Levy processes with variance-mean mixed normal distributions, in particular, 
to stable Levy processes, generalized hyperbolic and generalized variance-gamma Levy processes. 
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1 Introduction 

In financial mathematics the evolution of (the logarithms of) stock prices and financial indexes on small time 
horizons is often modeled by random walks. The simplest example of such an approach is the Cox-Ross- 
Rubinstein model (see, e. g., [25]). At the same time most successful (adequate) models of the dynamics of 
(the logarithms of) financial indexes on large time horizons are subordinated Wiener processes (processes of 
Brownian motion with random time such as generalized hyperbolic processes, in particular, variance gamma 
(VG) processes and normal\\inverse Gaussian (NIG) processes. These models very well describe the observed 
heavy-tailedness and leptokurticity of the empirical (statistical) distributions of the increments of financial 
indices and, in particular, of stock prices on comparatively short time intervals. 

Functional limit theorems are a quite natural link between random walks and subordinated Wiener processes. 
The operation of subordination gives a good explanation of the presence of heavy tails in the empirical 
distributions of the increments of (the logarithms of) stock prices and financial indexes. 

In the book m and the papers DUES] it was proposed to model the evolution of non-homogeneous chaotic 
stochastic processes, in particular, of the dynamics of stock prices and financial indexes, by random walks 
generated by compound doubly stochastic Poisson processes (compound Cox pocesses). A doubly stochastic 
Poisson process (also called a Cox process ) is a stochastic point process of the form Ni(A(t)), where Ni(t), t > 0, 
is a homogeneous Poisson process with unit intensity and the stochastic process A(t), t > 0, is independent of 
Ni(t) and possesses the following properties: A(0) = 0, P(A(f) < oo) = 1 for any t > 0, the sample paths of A(£) 
do not decrease and are right-continuous. A compound Cox process is a random sum of independent identically 
distributed random variables in which the number of summands follows a Cox process. Similar continuous-time 
random walks were considered in mm- 

This approach, based on the universal principle of non-decrease of entropy in closed systems, was developed 
in Esi uni uni- In the paper m this approach was successfully applied to modeling the evolution of limit order 
books in the high-frequency financial trading systems. In the framework of this approach the principal idea is 
that the moments at which the state of the system under consideration changes form a chaotic point stochastic 
process on the time axis. Moreover, this point process turns out to be non-stationary (time-non-homogeneous) 
because the changes of the state of the limit order book are to a great extent subject to the influence of 
non-stationary information flows. As is known, most reasonable probabilistic models of non-stationary (time- 
non-homogeneous) chaotic point processes are doubly stochastic Poisson processes also called Cox processes 
(see, e. g., IH0). These processes are defined as Poisson processes with stochastic intensities. Pure Poisson 
processes can be regarded as best models of stationary (time-homogeneous) chaotic flows of events 0. Recall 
that the attractiveness of a Poisson process as a model of homogeneous discrete stochastic chaos is due to 
at least two reasons. First, Poisson processes are point processes characterized by that time intervals between 
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successive points are independent random variables with one and the same exponential distribution and, as is well 
known, the exponential distribution possesses the maximum differential entropy among all absolutely continuous 
distributions concentrated on the nonnegative half-line with finite expectations, whereas the entropy is a natural 
and convenient measure of uncertainty. Second, the points forming the Poisson process are uniformly distributed 
along the time axis in the sense that for any finite time interval [t i, ^ 2 ], ti < ^ 2 ; the conditional joint distribution 
of the points of the Poisson process which fall into the interval [t \, £ 2 ] under the condition that the number of such 
points is fixed and equals, say, n, coincides with the joint distribution of the order statistics constructed from an 
independent sample of size n from the uniform distribution on [t \, £ 2 ] whereas the uniform distribution possesses 
the maximum differential entropy among all absolutely continuous distributions concentrated on finite intervals 
and very well corresponds to the conventional impression of an absolutely unpredictable random variable (see, 
e. g., mm)- 

In | I24| some functional limit theorems were proved establishing convergence of random walks generated by 
compound Cox processes with jumps possessing finite variances to Levy processes with symmetric distributions 
including symmetric strictly stable Levy processes. In the paper m these results were extended to a non- 
symmetric case and applied to modeling the evolution of the order flow imbalance process, an integral 
characteristic of the behavior of the limit order book. 

The present paper presents a further development of the models and techniques proposed in our previous 
papers (24, 27j. In this paper we improve and generalize the results of the mentioned papers by relaxing the 
conditions and correcting some inaccuracies. Our presentation essentially relies on the techniques developed in 

[Mi¬ 
ddle paper is organized as follows. Section 2 contains some basic material on the Skorokhod space, stable 
distributions and Levy processes. In section 3 we prove general functional limit theorem establishing the 
conditions for convergence of compound Cox processes to Levy processes in the Skorokhod space in terms 
of the behavior of the cumulative intensities of Cox processes. For this purpose we extend the classical results 
presented, say, in m • In Section 4 we consider the conditions for the convergence of compound Cox processes 
with elementary jumps possessing finite variances to the Levy processes with variance-mean mixed normal 
one-dimensional distributions, that is, to subordinated Wiener processes. 

2 Skorokhod space. Levy processes 

Let D = D[ 0,1] be the space of real functions defined on [0,1], right-continuous and having finite left-side limits. 

Let F be the class of strictly increasing mappings of [0,1] onto itself. Let / be a non-decreasing function on 
[0,1] with /(0) = 0, /(1) = 1. Set ||/|| = sup s#t jlog [(/(t) - f(s))/(t - s)] |. If ||/|| < 00 , then the function / is 
continuous and strictly increasing and, hence, belongs to J-. 

Define the distance do(x,y) in the set D[ 0,1] as the greatest lower bound of the set of positive numbers e, 
for which T contains a function / such that ||/|| ^ e and sup t | x(t) — | ^ e. 

It can be shown that the space D[0, 1] is complete with respect to the distance do. The metric space 
T> = (.D[0,1], do) is called the Skorokhod space. Everywhere in what follows we will consider stochastic processes 
as D-valued random elements. 

Let X : Xi,X 2 ,... be D-valued random elements. Let Tx be a subset of [0,1] such that 0 £ Tx , 1 £ Tx and 
if 0 < t < 1, then t £ Tx if and only if P(X(i) ^ X(t —)) = 0. The following theorem establishing sufficient 
conditions for the weak convergence of stochastic processes in V (denoted below as => and assumed asn-> 00 ) 
is well-known. 

Theorem A. Let (A'„(fi),..., X„(ffc)) => (X(t \),..., X(tk)) for any natural k and t \,..., tk belonging to Tx- 
Let P(A(1) ^ X (1—)) = 0 and let there exist a non-decreasing continuous function F on [0,1], such that for 
any e > 0 

P(| X n (t) - X n (t!)\ 2 e, \X n {t 2 ) - X n {t)\ >e) ^ e- 2v [F(t 2 ) - F(h)\^ (1) 

for ti ^ t ^ t 2 and n ^ 1, where v ^ 0, 7 > 1/2. Then X n => X. 

The PROOF of Theorem A can be found, for example, in [7]. 

Everywhere in what follows the symbol = stands for the coincidence of distributions. 

By a Levy process we will mean a stochastic process X(t), t ^ 0, possessing the following properties: (i) 
X(0) = 0 almost surely; (ii) X(t) is the process with independent increments, that is, for any N ^ 1 and 
to,ti, ...,£jv (0 < to ^ < ... < tx) the random variables X(t, 0 ), X{ti) — A(t 0 ), X{tx) — X(tN-i) are jointly 

independent; (iii) X(t) is a homogeneous process, that is, X(t+h)—X(t) = X(s+h)—X(s) for any s, t, h > 0; (iv) 
the process X(t) is stochastically continuous, that is, for any t ^ 0 and e > 0 lim s ^ t P(|X(t) — X(s)| > e) = 0; 
(v) sample paths of the process X(t) are right-continuous and have finite left-side limits. 
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Denote the characteristic function of the random variable X(t) as ipt{s) (V’t(s) = Ee lsX ^\ s £ R). The 
following statement describes a well-known property of Levy processes. 

Lemma 1. Let X = X(t), t ^ 0, be a Levy process. For any t > 0 the characteristic function of the random 
variable X (t) is infinitely divisible and has the form 

Ms) = [Ms)] t =[ Ee isX ^]\ set. (2) 

Conversely, let Y be an arbitrary infinitely divisible random variable. Then the family of infinitely divisible 
distributions with characteristic functions of the form [Ee lsy ] completely determines finite-dimensional 
distributions of a Levy process X(t), t Js 0, moreover, X(l) = Y. 

The properties of Levy processes are described in detail in El ESI- 

By G a fi(x) we will denote the distribution function of the strictly stable law with the characteristic exponent 
a and parameter 9 corresponding to the characteristic function g Qi 6i(s) = exp { —|s|“ exp { —|7r#asigns}}, s £ R, 
where 0 < a ^ 2, \6\ ^ 9 a = min{l, ^ — 1}. To symmetric strictly stable distributions there corresponds the 
value 9 = 0. To one-sided stable distributions there correspond the values 9 = 1 and 0 < a ^ 1. 

In what follows any random variable with the distribution function G at g(x), 0 < a < 2, will be denoted 
Z a ^g. It is well known that E\Z a: g\ s < oo for any S € (0,a), but the moments of orders higher or equal to a of 
the random variable Z ai g do not exist (see, e. g., pH]). 

The distribution function of the standard normal law (a = 2, 9 = 0) will be denoted <f>(:r), <h(a;) = 

HooVi^dz, <p(x) = 7=e-* 2 / 2 . 

It is well known that the distribution function G a ,o(x) of the symmetric strictly stable law can be represented 
as a scale mixture of normal laws: 


G a , 0 (x) 



xGl 


(see, e.g., ED, Theorem 3.3.1). To this representation there corresponds the following relation in terms of 
characteristic functions: 


0«,o(s) 



2 

s z u 

~~ 2 ~~ 


} 


dG a /2,i(u), 


s€l. 


A Levy process X(t), t ^ 0, will be called a-stable, if P(X(1) < x) = G a p(x), x € M. It can be shown (see, 
e.g., [ 5 ]) that if X(t), t ^ 0, is a Levy process, then X(t) is a-stable if and only if 


X(t)=t 1/a X{ 1), 0. 


( 3 ) 


3 Convergence of compound Cox processes to Levy processes 

In what follows without noticeable loss of generality we will consider stochastic processes defined for 0 ^ t ^ 1. 
Actually, this means that we consider the behavior of compound Cox processes on finite time horizons. The 
equality of the right bound of the horizon to one can be achieved by an appropriate choice of the units of 
measurement of time. In other words, we will concentrate on studying the case of the Skorokhod space T>. 

In order to introduce reasonable asymptotics which formalizes the condition of «infinite» growth of 
intensities of the flows of informative events, and makes it possible to construct asymptotic («heavy-traffic») 
approximations to Cox processes, fix a time instant t and introduce an auxiliary parameter n. Everywhere in 
what follows the convergence will be meant as n —> oo unless otherwise specified. So, consider a sequence of 
compound Cox processes of the form 


Qn{t) = Y, 


N[ n \K n (t)) 


X n ,i, t ^ 0 , 


( 4 ) 


where {N[ n \t), t ^ 0} is a sequence of Poisson processes with unit intensities; for each n = 1,2,... the 
random variables X n ,i, X n p,,... are identically distributed; for any n ^ 1 the random variables X n> i,X n> 2 ,... 
and the process N[ n \t), t ^ 0, are independent; for each n = 1,2,... A n (t), t ^ 0, is a subordinator, that is, a 
non-decreasing positive Levy process, independent of the process 


Z r 


At )" 5 (i) 


(t) = Y \ o, 

L — j i—i 


( 5 ) 
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and such that A„(0) = 0 . 

Assume that there exist 5 £ (0,1], 5\ ^ \ and the constants C n £ (0, oo) providing for all t £ (0,1] the 
validity of the inequality 

EA*(i) < (( C n t) s K ( 6 ) 

For example, assume that A n (t) is a stable Levy process, that is, P(A n (l) < x) = G ay i(x) with some 
0 < a < 1. Then for any p £ (0, a) C (0,1] we have EA£(1) < oo and, moreover, in accordance with (3) 
we have EA£(i) = t p/,a EA£(1), that is, condition (6) holds for any S £ [a/2, a) with = S/a £ [|, 1 ) and 

C n = (EA s n (l)) 1/a . 

Assume that 

0 < = E\X nyl f < oo. (7) 

for some /3 £ (0,1]. 

Everywhere in what follows for definiteness we assume that Y^!= i = 0. 

From (4) and (5) it is easy to see that Q n (t) = Z n (A n (t)). Since for each n ^ 1 both Z n (t) and A n (t) are 
independent Levy processes, and, moreover, A n (t) is a subordinator, then the superposition Q n (t) = Z n (A n (t)) 
is also a Levy process (see, e. g., Theorem 3.1.1 in mi). Hence the following statement follows. 

Lemma 2 . For any 0 ^ t\ < £2 < 00 and any n ^ 1 we have Q n {t 2) — Q n (t\) = Q n (t 2 — ii). 

Lemma 3. Let Q n (t) be a compound Cox process (4) satisfying conditions ( 6 ) and (7). Then for any t £ [0,1] 
and any e > 0 we have P(|Qn(i)| e) ^ ■ (Cut) 61 . 

PROOF. Since one-dimensional distributions of the Cox process (4) are mixed Poisson, we have 


P(\Q n (t)\>e) = P £■ 




J 3 =1 


x„ 


>')= E[L„ P(JV. ( "’(A~(t)) = *)P(| Xnj\ ><) = 




^ e 


dP(A„(t) < A). 


( 8 ) 


The change of the order of summation and integration is possible due to the obvious uniform convergence of the 
series. Continue ( 8 ) by the sequential application of the Markov and Jensen inequalities with 5 £ (0,1] taking 
part in ( 6 ) and /3 £ (0,1] taking part in (7). As a result we obtain 


P(IQ n (*)l > <0 



E OO \ A^' I ^—A k 

e~ x — E V 

k =0 kl I 1 




n,j 




dP(A n (t) < A) < 




X„ 



dP(A n (t) < A), 


(9) 


since with S £ (0,1] the function f(x) = x s is concave for x ^ 0. It is easy to see that E| ^2j=iXn,j\^ ^ 

Y2 k j=i E\X n ,j\P = kmf/ for 0 < /3 ^ 1. Therefore, continuing (9) with the account of the Jensen inequality for 
concave functions and ( 6 ), we obtain 


P( \Qn(f )| ^ e) ^ 


e ps 


( 00 




=0 


^)dP(A„(t) < A) = ^^^°°E[JV 1 ( 1 ) (A)] 4 dP(A n (t) < A) 


€ 


«) s 

f ps 


€ 


r oo ( / 3\<5 r°o 

/ [ElV 1 ( 1 ) (A)] 5 dP(A„(t) < A) = / X s dP(A n (t) < A) 

Jo e Jo 



• EA£(t) < 



• ( C n t) Sl . 


The lemma is proved. 

To establish weak convergence of the stochastic processes Q n (t) in the Skorokhod space 2?, first it is required 

to find the limit distribution of the random variables Q n (t) for each t > 0. The symbol —will denote 
convergence in distribution, that is, pointwise convergence of the distribution functions in all continuity points 
of the limit distribution function. 

Let t = 1. Denote N n = N[ n \ A„(l)). Assume that for some the convergence 


P{X nA + ... + X Uykn <x)^ H{x) 
takes place, where H ( x ) is some infinitely divisible distribution function. 


( 10 ) 
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Also assume that 


P(A„( 1 ) < k n x ) P (U < x), (11) 

where U is a nonnegative random variable such that its distribution is not degenerate in zero. Notice that since 
A n (t) is a Levy process, then the random variable U is infinitely divisible being the weak limit of infinitely 
divisible random variables. 

Lemma 4. Let N n = N^ n \ A„), n ^ 1, where {N^^, t 0}, n = 1,2,... are standard Poisson processes 
and A ra , n = 1,2,... are positive random variables such that for each n ^ 1 the random variable A n is 
independent of the process N[ n \t). Then P(N n < k n x) —A A(x) for some infinitely increasing sequence k n 
of real numbers and some distribution function A(x) if and only if P(A„ < k n x) A(x). 

For the proof see EDI- 

From Lemma 4 it follows that convergence (11) is equivalent to 


P(iVn < k n x) A P(J7 < x). 

By the Gnedenko-Fahim transfer theorem [91 conditions (10) and (12) imply that 


( 12 ) 


Qn(l)=X nA + ... + X n , Nn A Q, (13) 

where Q is a random variable with the characteristic function f(s) = / 0 °° (h(s)) u dP(U < w), h(s ) being the 
characteristic function corresponding to the distribution function H(x). Note that the distribution function 
H(x) may not satisfy the condition H(—x) = 1 — H(x) for all x ^ 0, that is, it may not be symmetric. 

Let Y be an infinitely divisible random variable with the distribution function H(x). Since both Y and U 
are infinitely divisible, we can define independent Levy processes Y(t) and 17(f), t 0, such that Y(l) = Y 
and 17(1) = Y. Then with the account of Lemma 1 it is easy to verify that f(s) = Ee ls ® = Eexp {isF(17(l))}, 

s £ R, that is, Q = Y(U(1)). Moreover, repeating the reasoning from [TB] (see Theorem 3.3.1 there), we can 
easily see that the random variable Q is infinitely divisible and hence, we can define a Levy process Q(t), t ^ 0, 
such that <3(1) = Q. From Lemma 1 and the abovesaid it follows that we can regard Q(t) as the superposition: 
Q(t)±Y(U{t)). 

Since according to (13) we have Q n ( 1 ) = J2i =i X n ,* ==> <3(1), and both Q n (t) and Q{t) are Levy processes, 
then, using (2) we can conclude that for any t > 0 


E JV»,i(A n (t)) j 

X n i — Q(t). 

1 


(14) 


Since the processes Q„(t) and Q(t), 0 ^ t ^ 1, are Levy processes, then almost all their sample paths belong 
to the Skorokhod space V. 

Consider the question what additional conditions are required to provide the weak convergence of the 
compound Cox process Q n (t ) to the Levy process Q(t) in the space T>. We will consider each of the conditions 
of Theorem A one by one. 

First, without loss of generality, let 0 ^ t± < t2 < ^ 1. The convergence (Q n {ti ),..., Q n {tk)) —A 

..., Q(tk )) is equivalent to the convergence 


{Qn(t l), Q n (t 2) 


Qnifk) Qnifk— l)) 


—*■ (Q(h),Q(t2) - ...,Q(t k ) - Q{t k - 1 )), (15) 

since the linear transform (x\,X 2 , ...,Xk-i,Xk) 1 —> (x±,X 2 — x±, ...,Xk — x k -i) of to is one-to-one and 
continuous in both directions. But convergence (15) follows from (14) and the fact that both Q n (t) and Q(t ) 
are Levy processes. 

Second, we have to check the condition P(Q(1) ^ Q(l—)) = 0. This condition holds if and only if 
lim t _»i_ P(|Q(1) — Q(t) | > e) = 0 for any e > 0 (see relation (15.16) in El)- Consider P(|Q(1) — Q(t)\ > e). 
Since Q(t) is a Levy process, then Q(l) — Q(t) = Q( 1 — t) by Lemma 2. Therefore, P(|Q(1) — Q(t)\ > e) = 
P(|Q(1 — t)\ > e). For each e > 0 and each t £ [0,1] there exists an e t & [e/2, e] such that the points ±e t are 

continuity points of the distribution function of the random variable Q( 1 — t). Since Q n (t) — Q(t) for each 
t £ [0,1], then P(|Q(1 — t)\ > e t ) = lim^oc P(|Q„(1 — t)| > e t ). Thus, for any e > 0 and any t £ [0,1] we have 

P(|Q(1 — i)| > e) ^ P(|Q(1 - t)\ > e t ) = lim P(|Q n (l - t)\ > e t ) < limsup P(|Q„(1 - t)\ > e t ) ■ (16) 
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Continuing (16) with the account of (6) and applying Lemma 3, for 6 £ (0,1] taking part in (6) we obtain 

P(|Q(1 - t) | > e) < limsup P(|<2„(1 - t)\ > e t ) < 

n—> oo 

< limsup (et P m^) S (C n \l -t |) 51 ((2/e) /3<5 |l - f|) dl limsup (ruffiC^ 1 . (17) 

n—> oo n—too 

Therefore, if 

K = limsupC^ 1 //| 5 ra^ < oo, (18) 

n—> OO 

then (17) implies 

lim P(|Q(1) - Q(t) | > e) < A(Ke~^) s lim_ |1 - t\ Sl = 0. 

Third, check condition (1) under the assumption that ( 6 ) and (18) hold. As it has been noted above, Q n (t) is 
a Levy process and hence, it has independent increments. Therefore, 

P(| Q n (t) - Q n {ti)\ > e, |Q„(£ 2 ) - Quit) | >e) = P(|Q„(t) - Q n (ti)\ ^ e) • P(|Q n (* 2 ) - Quit) | > e). (19) 

Consider the first multiplier on the right-hand side of (19). By Lemma 2, Q n (t) — Q n (t\) = Q„(t — ti). With 
the account of (18), by Lemma 3 we obtain 

P(| Q n (t) - Q n {t i)| >e)= P(| Quit - h) | > c) < (Ke~P) s \t - h\ Sl . (20) 

For the second multiplier on the right-hand side of (19) we similarly obtain 

P(|Qn(* 2 ) - Qn(t )| > e) = P(|Q„(i 2 - t )| > c) < (A"e _/3 )' 5 |t 2 - t\ s \ (21) 

Thus, from (20) and (21) it follows that 

P(|Qn(f) - Qn(ti)\ > e, I Q n (t 2 ) - Q n (f) \ > e) < (J^e-^) 25 [(t - h)(t 2 - t )]* 1 (22) 

It is easy to see that for any ^ t ^ t 2 we have (t — ti)(f 2 — f) ^ \(t 2 — t \) 2 . Substituting this estimate in (22) 

we obtain P(|<3n(£) — Qn(ti)\ ^ e, |Qn(i 2 ) — Q n (i)| > e) ^ e~ 2 ^ 6 [\K(t 2 — £i)] 2<Sl . Therefore, if conditions ( 6 ) 
and (18) hold, then condition (1) holds with F(t) = \Kt, v = /35 and 7 = Si. 

Summarizing this reasoning related to checking the conditions of Theorem A, we arrive at the following 
statement. 

Theorem 1. Assume that the random variables n = 1,2,..., (the jumps of the compound Cox 

process Q n (t), see (4)), satisfy conditions (10) with some k n £ N and (7) with some £ (0,1]. Let the processes 
Q n (t) be lead by non-decreasing positive Levy processes A n (t) satisfying conditions ( 6 ) with some S £ (0,1], 
<5i > i and (11) with the same k n . Also assume that condition (18) holds. Then the processes Q n (t) weakly 
converge in the Skorokhod space T> to the Levy process Q(t) such that 

fOO 

Eexp{*sQ(l)} = / (h(s)) u dP(U <u), set, (23) 

Jo 

where h(s) is the characteristic function corresponding to the distribution function H(x) in (10). 

It is worth noting that actually Theorem 1 deals with the well-studied weak convergence of special 
semimartingales with stationary increments, see, e. g., m- However, the superposition-type structure of the 
processes considered in the present paper makes it possible to relax the conditions required in the general case, 
say, in Corollary VII.3 .6 of [T5] where it is assumed that (in our terminology) <5 = <$1 = 1. Moreover, in [23] and 
some subsequent papers a more restrictive condition was used instead of (18). 

From Theorem 1 it obviously follows that if in (10) H(x) = G a p(x) for some admissible a £ (0,2) and 
9 £ [— 1 ,1] and the processes A n (t) are asymptotically degenerate, that is, in ( 11 ) P(U = 1) = 1, then in (23) 
Eexp{fsQ(l)} = fl ai e(s), that is, the limiting process is the stable Levy process. 

However, in m it was demonstrated that stable Levy processes can appear as limits for compound Cox 
processes even when the variances of elementary increments of a compound Cox process are finite. To generalize 
this result, consider the following corollary of Theorem 1. 

COROLLARY 1. Assume that the random variables {X n j}j^ 1 , n = 1,2,..., satisfy condition (10) with some 
k n £ N and H(x) = G a p(x) for some a £ (0,2] so that condition (7) holds for any ft £ (0, a) f] (0,1]. Let 
the compound Cox processes Q n (t) be lead by non-decreasing positive Levy processes A n (t) satisfying conditions 
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(11) with the same k n and P(U < x) = G Q yi(x) for some a' £ (0,1] so that condition (6) holds with some 5 £ 
[of / 2, c/) and <5i = 5/a 1 £ [i, 1). Also assume that condition (18) holds for some /3 £ (0, a). Then the processes 
Q n (t) weakly converge in the Skorokhod space T> to the Levy process Q(t ) such that P(<5(1) < x) = G aa ',o(x), 
x £ R. Proof. To prove this result is suffices to notice that in the case under consideration f a ,o(s) = e _ l s l , 
s £ R, so that in (23) 

pOO poo 

Eexp{7sQ(l)}= / e"“ |sr dG a /,i(u)= / e - |s “ 1/ “ |Q dG Q /,i(u) = Eexp {isZ^oZ 1 ^}, s £ R, 

Jo Jo 

where the random variables Z a $ and Z a ',i are independent. But from Theorem 3.3.1 in [31] it follows that 
ZafiZ'J'l = Z aa tfi. The corollary is proved. 

In turn, the result from |24| mentioned above follows from Corollary 1, if a = 2. This case corresponds to 
the situation in which the variances of the summands (elementary jumps) are assumed finite. As it has been 
already said, in most applied problems there are no reasons to reject this assumption. Therefore in what follows 
we will concentrate our attention on the case where the elementary increments of the compound Cox process 
have finite variances and consider the conditions of convergence of compound Cox processes to some popular 
models, in particular, to Levy processes with variance-mean mixed normal one-dimensional distributions such 
as generalized hyperbolic Levy processes or generalized variance-gamma Levy processes. 

4 Levy processes with variance-mean mixed normal distributions as 
asymptotic approximations to compound Cox processes 

Denote a n = EX^p and a% = DX n i. From the classical theory of limit theorems it is well known that if, as 
n —y oo, the conditions 

knan t a, k n a n ^ o and k n P(X n ^i Un) I(|AVi,i a n \ jz- e) y 0 (24) 

hold for some a £ R, 0 < a 2 < oo and any e > 0, then convergence (10) takes place with H{x ) = $(cr _1 (a; — a)). 
In this case the distribution function F{x) of the limit random variable Q(l) in Theorem 1 is a variance-mean 
mixture of normal laws. Recently it was demonstrated that normal variance-mean mixtures appear as limiting in 
simple limit theorems for random sums of independent identically distributed random variables m ■ Namely, let 
n = 1, 2 ,..., be a double array of row-wise (for each fixed n) independent and identically distributed 
random variables. Let {v n }n^i be a sequence of integer nonnegative random variables such that for each n ^ 1 
the random variables u n , 2 , ■ ■ ■ are independent. Denote S n< k = £»,i + ... + £, n ,k- The following theorem 

was proved in EH- 

Theorem B. Assume that there exist: a sequence {k n } n ^ i of natural numbers and finite numbers a £ R 
and a > 0 such that 

P(S n , kn <x) (25) 

Assume that v n —^ oo in probability. Then the distribution functions of random sums S n ^ n converge to some 

distribution function F(x) : P(S’ n ,j/„ < ——* F(x), if and only if there exists a distribution function A(x) such 

that 

A( 0) = 0, F(x) = $(JL^)dA(z), 

and P (y n < xk n ) A(x). 

Theorem B and Lemma 3 yield the following result. 

Theorem 2. Assume that the random variables {X n j}j^i, n = 1,2,..., (the jumps of the compound Cox 
process Q n (t), see (4)) possess finite variances and satisfy conditions (24) with some k n £ N, a £ M and a > 0. 
Let the the processes Q n {t) be lead by non-decreasing positive Levy processes A n (t) satisfying condition (6) with 
some S £ (0, 1], 6i ^ i. Also assume that 

K = lira sup Gf l/<5 (cr„ + |a n |) < oo. (26) 

n—> oo 

Then the processes Q n (t) weakly converge in the Skorokhod space T> to a Levy process Q(t) if and only if there 
exists a nonnegative random variable U such that 

P(Q(l)<x)= [ $(?—^)dP{U<u), £ £ R, (27) 

J o v <?V U ' 
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and condition (11) holds with the same k n . 

PROOF. It suffices to take /3 = 1 in the proof of lemma 3 and notice that m* ^ cr„ + \a n \ so that condition 
(26) can play the role of (18). 

The class of distributions of form (27) was systematically considered by O. Barndorff-Nielsen and his 
colleagues HHIl] in order to introduce generalized hyperbolic distributions and study their properties. 

The class of normal variance-mean mixtures (27) is very wide. For example, it contains generalized 
hyperbolic laws with generalized inverse Gaussian mixing distributions mm, generalized variance gamma (GVG) 
distributions with generalized gamma mixing distributions |221123] , symmetric strictly stable laws with strictly 
stable mixing distributions concentrated on the positive half-line, generalized exponential power distributions, 
and many other types. 

Generalized hyperbolic distributions demonstrate exceptionally high adequacy when they are used to describe 
statistical regularities in the behavior of characteristics of various complex open systems, in particular, turbulent 
systems and financial markets. There are dozens of dozens of publications dealing with models based on 
generalized hyperbolic distributions. Recently it was discovered that generalized variance gamma distributions 
demonstrate even better fit to empirical data. Therefore below we will concentrate our attention on functional 
limit theorems establishing the convergence of compound Cox processes to generalized hyperbolic Levy processes 
and generalized variance gamma Levy processes yielding the possibility of the use of such processes as convenient 
«heavy-traffic» asymptotic approximations. 

Denote the density of the generalized inverse Gaussian distribution by pgig(x\ C hi A), 


Pgig(x;v, p, A) 


A ^/ 2 

2pyGK v {y[p\) 


x -*. ex p{_i(£ + A*)}, 


X > 0. 


Here p > 0, A ^ 0 if < 0; /x > 0, A > 0 if z/ = 0 and /x ^ 0, A > 0 if v > 0, K„{z) is the modified Bessel function 
of the third kind with index i/, K v (z) = \ / 0 °° y U ~ 1 exp { — § {y + }dy, z € C, Re z > 0. The corresponding 
distribution function will be denoted Pgig{x', v, hi A). 

The class of generalized inverse Gaussian distributions is rather rich and contains, in particular, both 
distributions with exponentially decreasing tails (gamma-distribution (/i = 0, v > 0)), and distributions 
whose tails demonstrate power-type behavior (inverse gamma-distribution (A = 0, v < 0), inverse Gaussian 
distribution [v = — L) and its limit case as A —► 0, the Levy distribution (the stable distribution G i x (x) with 
the characteristic exponent equal to i and concentrated on the nonnegative half-line, the distribution of the 
time for the standard Wiener process to hit the unit level)). 

In 1977-78 O. Barndorff-Nielsen [1)[5] introduced the class of generalized hyperbolic distributions as the class 
of special normal variance-mean mixtures. For convenience, we will use a somewhat simpler parametrization. 
Let a £ R, a > 0. If the generalized hyperbolic distribution function with parameters a, a, zx, p, A is denoted 
Pgh{x] a, cr, z/, /x, A), then by definition, 


I ( x — olz \ 

P G H(x-,a,cr,i>,h, A) = / $ - -yf- )pgig(z',Vi hi A)ck, x € R. (28) 

J 0 V oyjz / 

Note that in (28) mixing is carried out simultaneously by both location and scale parameters, but since these 
parameters are directly linked in (28), then actually (28) is a one-parameter mixture. 

It is well-known that all generalized hyperbolic distributions are infinitely divisible mia. 

From Theorem 2, Corollary 1 and Theorem B with the account of the equivalence of relations (11) and (12) 
we easily obtain the following result on the convergence of compound Cox processes to generalized hyperbolic 
Levy processes. 

COROLLARY 2. Assume that the random variables {Xn,j}j^l, n = 1, 2,..., ( the jumps of the compound Cox 
process Q n (t), see (4)) possess finite variances and satisfy conditions (24) with some k n £ N, a G K and a > 0. 
Let the the processes Q n {t) be lead by non-decreasing positive Levy processes A n (t) satisfying condition (6) with 
some S 6 (0,1], ^ Also assume that condition (26) holds. Then the processes Q n (t) weakly converge in the 

Skorokhod space T> to a generalized hyperbolic Levy process Q(t) such that P(Q(1) < x) = PGH(x',a,a,u, y, X) 

if and only if P(A„(1) < k n x) — Pgig{x\ v , y,, A) with the same k n , v, p, and A. 

The class of generalized gamma (GG) distributions was introduced in the paper |30| . Any representative of 
this class is defined by the probability density 


6 k "T(K) 


x exp 


{-(in. 


Pgg{x;v,k,5) 


x > 0 





where the parameters v £ R, k, S £ R + respectively determine the power, shape and scale of the generalized 
gamma distribution. Here T(k) is Euler’s gamma-function. 

The family of GG-distributions contains practically all most popular absolutely continuous distributions 
concentrated on R + . In particular, this class contains (i) the gamma-distribution (v = 1), (ii) the exponential 
distribution (v = 1, k = 1), (iii) the Erlang distribution ( v = 1, n £ N), (iv) the chi-square distribution 
{y = 1, 6 = 2), (v) the Nakagami distribution (v = 2), (vi) the half-normal distribution (the distribution of 
the absolute value of the standard normal random variable) (v = 2, k = ^), (vii) the Rayleigh distribution 
(y = 2, k = 1), (viii) the chi-distribution (y = 2, 8 = a/2) , (ix) the Maxwell distribution (y = 2, k = 3/2), 
(x) the Weibull-Gnedenko distribution (k = 1), (xi) the inverse gamma-distribution (v = —1), (xii) the Levy 
distribution, (xiii) the lognormal distribution (k —> oo). 

Unlike generalized inverse Gaussian laws, generalized gamma laws contain the distributions with exponential 
power type of decrease of the tail where the posiltive exponent power may be arbitrary. 

In (22j the family of generalized variance gamma (GVG) distributions was introduced as the family of normal 
variance-mean mixtures of the form 


Pgvg{x; a, a, v,k,5) 



$ 


x — au 

(Ty/U 


PGG{u\v,n,S)du, 


x £ R, 


Both classes, GH and GVG distributions, contain (a) symmetric and non-symmetric (skew) Student 
distributions (including Cauchy distribution), to which in (32) there correspond inverse gamma mixing 
distributions; (6) variance gamma (VG) distributions) (including symmetric and non-symmetric Laplace 
distributions), to which in (32) there correspond gamma mixing distributions; (c) normal\\inverse Gaussian 
(NIG) distributions to which in (32) there correspond inverse Gaussian mixing distributions. However, the class 
of GVG laws includes normal mixtures with tails decreasing as e - ^ with 0 < v < 1, which are not included 
in GH-distributions. This type of tail behavior is of great practical interest. 

Not all of GVG distributions are infinitely divisible. But as concerns the important case of Weibull-type 
decreasing tails mentioned above, it turns out that for v £ (0,1] these laws are infinitely divisible. Recall that in 
the notation introduced above the Weibull-Gnedenko distribution corresponds to the GG-density pgg{%\ v, 1,5). 

Lemma 5. If v ^ 1, then the Weibull-Gnedenko distribution is infinitely divisible. 

PROOF. As was shown in the paper [26], the Weibull-Gnedenko distribution with u ^ 1 is mixed exponential. 
But in m it was proved that all mixed exponential laws are infinitely divisible. The lemma is proved. 

According to lemma 5, if v ^ 1, then a Levy process 17(f), t ^ 0, can be defined so that P(t/(1) < x) = 
Pgg{x ; v , 1, 6), itl. Such a process will be called the Levy-Weibull process. Let W(t) be a standard Wiener 
process independent of U(t). 

COROLLARY 3. Assume that the random variables {X n ,j}j^i, n = 1,2, ..., (the jumps of the compound Cox 
process Q n (t), see (4)) possess finite variances and satisfy conditions (24) with some k n £ N, a = 0 and a > 0. 
Let the the processes Q n (t) be lead by non-decreasing positive Levy processes A n (t) satisfying condition (6) with 
some 5 £ (0,1], Also assume that condition (26) holds. Then the processes Q n (t ) weakly converge in the 

Skorokhod space V to a subordinated Wiener process W(U(t)) with the subordinator U(t) being the Levy-Weibull 

process with u ^ 1 if and only if P(A„(1) < k n x) — Pqg{x\ v , 1, S) with the same k n . 
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